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Axions solve the Strong CP Problem and are a cold dark matter candidate. The combined 
constraints from accelerator searches, stellar evolution limits and cosmology suggest that the 
axion mass is in the range 3 • 10"^ > nia > 10'^' eV. The lower bound can, however, be 
relaxed in a number of ways. I discuss the constraint on axion models from the absence of 
isocurvature perturbations. Dark matter axions can be searched for on Earth by stimulating 
their conversion to microwave photons in an electromagnetic cavity permeated by a magnetic 
field. Using this technique, limits on the local halo density have been obtained by the Axion 
Dark Matter eXperiment. 



1 Introduction 

The standard model of elementary particles has been an extraordinary breakthrough in our 
description of the physical world. It agrees with experiment and observation disconcertingly 
well and surely provides the foundation for all further progress in our field. It does however 
present us with a puzzle. 

Indeed the action density includes, in general, a term 

■^stand mod = ••• "I" '^2tt^^ ^'^^ ^ 

where G^^, are the QCD field strengths, g is the QCD coupling constant and ^ is a parameter. 
The dots represent all the other terms in the action density, i.e. the terms that lead to the 
numerous successes of the standard model. Eq. perversely shows the one term that isn't 
a success. Using the statement of the chiral anomaly^, one can show three things about that 
term. First, that QCD physics depends on the value of the parameter 6 because in the absence 
of such dependence QCD would have a symmetry in the chiral limit, and we know QCD 



has no such ^^(1) symmetry^. Second, that 9 is cycHc, that is to say that physics at 9 is 
indistinguishable from physics at ^ + 2it. Third, that an overall phase in the quark mass matrix 
ruq can be removed by a redefinition of the quark fields only if, at the same time, 9 is shifted to 
9 — arg det rUq. The combination of standard model parameters = — arg det rUq is independent 
of quark field redefinitions. Physics therefore depends on solely through 0. 

Since physics depends on 9, the value of 9 is determined by experiment. The term shown 
in Eq. violates P and CP. This source of P and CP violation is incompatible with the 
experimental upper bound on the neutron electic dipole moment unless \9\ < 10^^'^. A new 
improved upper limit on the neutron electric dipole moment i\dn\ < 3.0 • 10^^^ e cm) was 
reported at this conference by P. Geltenbort^^. The puzzle aforementioned is why the value of 
9 is so small. It is usually referred to as the "Strong CP Problem". If there were only strong 
interactions, a zero value of 9 could simply be a consequence of P and CP conservation. That 
would not be much of a puzzle. But there are also weak interactions and they, and therefore 
the standard model as a whole, violate P and CP. So these symmetries can not be invoked to 
set ^ = 0. More pointedly, P and CP violation are introduced in the standard model by letting 
the elements of the quark mass matrix niq be arbitrary complex numbers^. In that case, one 
expects arg det mq, and hence 9, to be a random angle. 

The puzzle is removed if the action density is instead 



1 (/^ (2(x) 

•^stand mod + axion — • • • 



where a(x) is a new scalar field, and the dots represent the other terms of the standard model, fa 
is a constant with dimension of energy. The aG ■ G interaction in Eq. ^ is not renormalizable. 
However, there is a recipe for constructing renormalizable theories whose low energy effective 
action density is of the form of Eq. ((2}. The recipe is as follows: construct the theory in such a 
way that it has aU{l) symmetry which (1) is a global symmetry of the classical action density, 
(2) is broken by the color anomaly, and (3) is spontaneously broken. Such a symmetry is called 
Peccei-Quinn symmetry after its inventors'^. Weinberg and Wilczek- pointed out that a theory 
with a [/pq(1) symmetry has a light pseudo-scalar particle, called the axion. The axion field is 
a{x). fa is of order the expectation value that breaks C/pq(1), and is called the "axion decay 
constant" . 

In the theory defined by Eq. 9 = ~ '^^^ ^.rg depends on the expectation value of 
a{x). That expectation value minimizes the effective potential. The Strong CP Problem is then 
solved because the minimum of the QCD effective potential V{9) occurs at 9 = 0^. The weak 



interactions induce a small value for 9^^, of order 10"^'^, but this is consistent with experiment. 

The notion of Peccei-Quinn (PQ) symmetry may seem contrived. Why should there be a 
U{1) symmetry which is broken at the quantum level but which is exact at the classical level? 
However, the reason for PQ symmetry may be deeper than we know at present. String theory 
contains many examples of symmetries whi ch are ex act classically but which are broken by 
quantum anomalies, including PQ symmetry' ' -^ l -^^ l. Within field theory, there are examples 
of theories with automatic PQ symmetry, i.e. where PQ symmetry is a consequence of just the 
particle content of the theory without adjustment of parameters to special values. 

The first axion models had fa of order the weak interaction scale and it was thought that 
this was an unavoidable property of axion models. However, it was soon pointed out that 
the value of fa is really arbitrary, that it is possible to construct axion models with any value 
fa- A value of fa far from any previously known scale need not lead to a hierarchy problem 
because PQ symmetry can be broken by the condensates of a new technicolor-like interaction 

The properties of the axion can be derived using the methods of current algebra^l. The 



axion mass is given in terms of fa by 

lO^GeV 

ma-6eV . (3) 

J a 

All the axion couplings are inversely proportional to fa- For example, the axion coupling to two 
photons is: 

a a(x) ^ 

Cayy = "^7 T^E -B . (4) 

J a 

Here E and B are the electric and magnetic fields, a is the fine structu re c onstant, and is a 
model-dependent coefficient of order one. = 0.36 in the DFSZ model^^l whereas g-y = -0.97 
in the KSVZ model The coupling of the axion to a spin 1/2 fermion / has the form: 



^ajf = ^df^aflbf (5) 



fa 

where gj is a model-dependent coefficient of order one. In the KSVZ model the coupling to 
electrons is zero at tree level. Models with this property are calle d 'h adronic'. 

The axion has been searched for in many places but not found^l. The resulting constraints 
may be summarized as follows. Axion masses larger than about 50 keV are ruled out by particle 
physics experiments (beam dumps and rare decays) and nuclear physics experiments. The next 
range of axion masses, in decreasing order, is ruled out by ste llar evolution arguments. The 
longevity of red giants rules out 200 keV > rria > 0.5 in the case of hadronic axions, 

and 200 keV > nia > 10~^ eV^^^ in the case of axions with a large coupling to electrons 
[ge = 0(1) in Eq. The duration of the neutrino pulse from Supernova 1987a rules out 2 
eV > nia > 3 • 10~^ eV^^. Finally, there is a lower limit, ^ 10~^ eV, from cosmology 
which will be discussed in detail in the next section. This leaves open an "axion window": 
3 • 10^"^ > nia ^ 10^^ eV. We will see, however, that the lower edge of this window (10^^ eV) 
is much softer than its upper edge. 



2 Axion cosmology 

The implications of the existence of an axion for the history of the early universe may be briefly 
described as follows. At a temperature of order /„, a phase transition occurs in which the 
Upq{1) symmetry becomes spontaneously broken. This is called the PQ phase transition. At 
these temperatures^the non-perturbative QCD effects which produce the effective potential 
V{6) are negligible the axion is massless and all values of {a{x)) are equally likely. Axion 
strings appear as topological defects. One must distinguish two scenarios, depending on wether 
inflation occurs with reheat temperature lower (case 1) or higher (case 2) than the PQ transition 
temperature. In case 1 the axion field gets homogenized by inflation and the axion strings are 
'blown' away. 

When the temperature approaches the QCD scale, the potential V{0) turns on and the axion 
acquires mass. There is a critical time, defined by ma{ti)ti = 1, when the axion field starts to 
oscillate in response to the turn-on of the axion mass. The corresponding temperature Ti ~ 1 
GeV'231. The initial amplitude of this oscillation is how far from zero the axion field lies when 
the axion mass turns on. The axion field oscillations do not dissipate into other forms of energy 
and hence contribute to the cosmological energy density today . This contribution is called 
of 'vacuum realignment'. It is further described below. Note that the vacuum realignment 
contribution may be accidentally suppressed in case 1 if the homogenized axion field happens 
to lie close to zero. 

In case 2 the axion strings radiate axion J^^^ from the time of the PQ transition till ti when 
the axion mass turns on. At ti each string becomes the boundary of domain walls. If = 1, 



the network of walls bounded by strings is unstable* * ' and decays away. If > 1 there is 
a domain wall problem"^ because axion domain walls end up dominating the energy density, 
resulting in a universe very different from the one observed today. There is a way to avoid this 
problem by introducing an interaction which slightly lowers one of the N vacua with respect to 
the others. In that case, the lowest vacuum takes over after some time and the domain walls 
disappear. There is little room in parameter space for that to happen and we will not consider 
this possibility"^ further here. Henceforth, we assume = 1. 

In case 2 there are thre e contributions to the axion cosmological energy density. One contri- 
butionl^ '^ l ^^ l ^^ l ^-*- ! ^^ ! ^^ ! ^ -^ is from axions that were radiated by axion strings before ti. A second 
contribution is from axions that were produced in the decay of walls bounded by strings after 
35 36 29 ^ third contribution is from vacuum realignment 

Let me briefly indicate how the vacuum alignment contribution is evaluated. Before time ti, 
the axion field did not oscillate even once. Soon after ti, the axion mass is assumed to change 
sufficiently slowly that the total number of axions in the oscillations of the axion field is an 
adiabatic invariant. The average number density of axions at time ti is 

Mh) ^ lma{ti){a^{ti)) ~ TT/f i (6) 

In Eq. ©, we used the fact that the axion field a{x) is approximately homogeneous on the 
horizon scale ti. Wiggles in a{x) which entered the horizon long before ti have been red-shifted 
awayl^. We also used the fact that the initial departure of a{x) from the nearest minimum is 
of order fa- The axions of Eq. Q are decoupled and non-relativistic. Assuming that the ratio 
of the axion number density to the entropy density is constant from time t^ till today, one finds 

^^'^ - 2 Uoi2GeV j [hj 
for the ratio of the axion energy density to the critical density for closing the universe, h is the 
present Hubble rate in units of 100 km/s.Mpc. The requirement that axions do not overdose 
the universe implies the constraint nia ^ 6 • 10~^ eV. 

The contribution from axion string decay has been debated over the years. The main issue 
is the energy spectrum of axions radiated by axion strings. Battye and Shellard'"^^' have carried 
out computer simulations of bent strings (i.e. of wiggles on otherwise straight strings) and have 
concluded that the contribution from string decay is approximately ten times larger than that 
from vacuum realignment, implying a bound on the axion mass approximately ten times more 
severe, say rua ^ 6 • 10~^ eV instead of rua ^ 6 • 1 0~^ eV . My collaborators and I have done 
simulations of bent strings^^, of circular string loops^^^^and non-circular string loops We 
conclude that the string decay contribution is of the same ord er of magnitude than that from 
vacuum realignment. Yamaguchi, Kawasaki and Yokoyama^^ have done computer simulations 
of a network of strings in an expanding universe, and concluded that the contribution from 
string decay is approximately three times th at o f vacuum realignment. The contribution from 
wall decay has been discussed in detail in ref.l^. It is probably subdominant compared to the 
vacuum realignment and string decay constributions. 

It should be emphasized that there are many sources of uncertainty in the cosmological 
axion energy density aside from the uncertainty about the contribution from string decay. The 
axion energy density may be diluted by the entropv release from heavy particles which decouple 
before the QCD epoch but decay afterwards!^, or by the entropy release associated with a 
firs t order QCD phase transition. On the other hand, if the QCD phase transition is first order 

an abrupt change of the axion mass at the transition may increase fi^. If inflation occurs 
with reheat temperature less than Tpg, there may be an accidental suppression of 0,a because 
the homogenized axion field happens to lie close to a CP conserving minimum. Because the 



RHS of Eq. (7) is multiplied in this case by a factor of order the square of the initial vacuum 
misalignment angle ^^^j^^ which is randomly chosen between — vr and +7r, the probability that 



fla is suppressed by a factor x is of order y/x. Recently, Kaplan and Zurek proposed a model 
in which the axion decay constant fa is time-dependent, the value fa{ti) during the QCD 
phase-transition being much smaller than the value fa today. This yields a suppression of the 
axion cosmological energy density by a factor ( -^"j*^^ )^ compared to the usual case [replace fa 
by fa{ti) in Eq. ©]. Finally, the axion density may be diluted by 'coherent deexcitation', i.e. 
adiabatic level crossing of rnait) with the mass of some other pseudo-Nambu-Goldstone boson 
which mixes with the axion 

The axions produced when the axion mass turns on during the QCD phase transition are cold 
dark matter (CDM) because they are non-relativistic from the moment of their first appearance 
at 1 GeV temperature. Studies of large scale structure formation support the view that the 
dominant fraction of dark matter is CDM. Any form of CDM necessarily contributes to galactic 
halos by falling into the gravitational wells of galaxies. Hence, there is excellent motivation to 
look for axions as constituent particles of our galactic halo. 

There is a particular kind of dumpiness which affects axion dark matter if there is no 
inflation after the Peccei-Quinn phase transition (case 2). This is due to the fact that the dark 
matter axions are inhomogeneous with 5p/ p ~ 1 over the horizon scale at temperature Ti ~ 1 
GeV, when they are produced at the start of the QCD phase-transition, combined with the fact 
that their velocities are so small that they do not erase these inhomogeneities by free-streaming 
before the time t^q of equality between the matter and radiation energy densities when matter 
perturbations can start to grow. These particular inhomogeneities in the axion dark matter 
are in the n on-li near regime immediately after time tgg and thus form clumps, called 'axion 
mini-clusters '1^. They have mass M^c — 1Q~^^Mq and size Imc — 10^^ cm. 

3 Axion isocurvature perturbations 

If inflation occurs after the Peccei-Quinn phase transition, i.e. if the reheat temperature after 
inflation Trh is less than the temperature Tpq at which f7pQ(l) is restored (case 1), the quantum 
mechanical fl uctuati ons of the axion field during the inflationary epoch cause isocurvature density 
perturbations l ^*^ ! ^^ ! in the early universe. The cosmic microwave background obervations are 
consistent with purely adiabatic density perturbations and therefore place a constraint, which 
we now discuss. 

Fluctuations generated during inflation in a massless weakl y co upled scalar field, such as the 
inflaton or the axion, are characterized by the power spectrum!^ 



where x are comoving coordinates. Eq. is often written in the shorthand notation 5a = 
The axion field fluctuations are "frozen" from the time their wavelengths exceed the horizon size 
HJ^ during the inflationary epoch till the time their wavelenghs reenter the horizon long after 
inflation has ended. 

At the start of the QCD phase transition, the local value of the axion field a(x, t) determines 
the local number density of cold axions produced by the vacuum realignment mechanism 



where a(x, ti) = a{x,ti)/fa is the local misalignment angle. The fluctuations in the axion field 






(9) 



produce perturbations in the axion dark matter density 



where ai = a{ti) = faCti is the initial value of the axion field, at the start of the QCD phase 
transition, common to our entire visible universe. These perturbations initially obey dp^" = 
—5pf" since the vacuum realignment mechanism converts energy stored in the quark-gluon 
plasma into axion rest mass energy. Such perturbations are commonly called "isocurvature 
perturbations" because they do not initially produce a source for the Newtonian potential $ = 
^(500 ~ !)• Note that in case 1 the density perturbations in the cold axion fluid have both 

adiabatic and isocurvature components. The adiabatic perturbations = = ^) are 

produced by the quantum mechanical fluctuations of the inflaton field during infiation, whereas 
the isocurvature perturbations [Spa"{ti) — —6pr{tiy^°] are produced by the quantum mechanical 
fiuctuations of the axion field during that same epoch. 

Isocurvature perturbations make a different imprint on the cosmic microwave background 
than do adiabatic ones. The CMBR obs ervations are consistent with pure adiabatic perturba- 
tions. According to P. Crotty et al. the fraction of cold dark matter perturbations which 
are isocurvature can not be larger than 31%. This places a constraint on axion models if the 
Peccei-Quinn phase transition occurs before infiation (case 1). Allowing for the possibility that 
only part of the cold dark matter is axions, the CMBR constraint of ref. 1^ is 

6p^- _ 6p^- Pa _ Hj na ^ Q 3^ Sprn ^ ^^^^ 



PCDM Pa PCDM T^faOil ^CDM Pn 

where we used Eq. ()1U() . is the amplitude of the primordial spectrum of matter perturbations. 
It is related to the amplitude of large scale (low multipole) CMBR anisotropies through the 
Sachs- Wolfe effecfS^. The observations imply'® ^ ~ 4.6 10"^ 

In terms of ai, the cold axion energy density is given in case 1 by 

where we assumed h ~ 0.7. It has been remarked by many authors, starting with S.-Y. Pi^^, 
that it is possible for fa to be much larger than lO^^GeV because ai may be accidentally small 
in our visible universe. The requirement that < r^cDM = 0.22"^ implies 

Since — vr < ai < +7r is the a-priori range of ai values and no particular value is preferred over 
any other, |^| may be taken to be the "probability" that the initial misalignment angle has 
magnitude less than |ai|. If | — | = 4 • 10~^, for example, fa may be as large as 10^^ GeV. 

The presence of isocurvatme perturbations constrains the small a, scenano m two waysEi. 
First, it makes it impossible to have ai arbitrarily small since 

ai > 6ai = . (14) 

2vr/a 

Combining Eqs. ((T^ and (fTl)) . we obtain the bound 

A/<e-10-Gev(^)* . (15) 



Second, one must require axion isocurvature perturbations to be consistent with CMBR obser- 
vations. Combining Eqs. (|TT|) and (|T^ . and setting JlcDM = 0.22, = 4.6 10~^, we obtain 

A,<10-GeV (^^)* . (16) 

Let us keep in mind that the bounds (|15() and ()16() pertain only if Trh < TpQ. One may, for 
example, have = 0.22, ~ lO^^ Qgy, and A/ ~ 10^^ GeV, provided Trh Z 10^^ GeV, which 
is possible if reheating is sufficiently efficient. 



4 Dark matter axion detection 



An electromagnetic cavity permeated by a strong static magnetic field can be used to detect 
galactic halo axions^^. The relevant coupling is given in Eq. Galactic halo axions have 
velocities /3 of order 10~^ and hence their energies Ea = ma + ^iriaP'^ have a spread of order 
10~^ above the axion mass. When the frequency uj = Inf of a cavity mode equals nia, galactic 
halo axions convert resonantly into quanta of excitation (photons) of th at cavity mode. The 
power from axion photon conversion on resonance is found to be^^^^l 



V BipaC—Mm{QL,Qa) 

faJ rUa 



0.5 10-26Watt ^ ^ \f^o 



,500 liter/ V7 Tesla 
Pa \ f m. 



10-''ci^/ V27r(GHz); 



]Min{QL,Qa) (17) 



where V is the volume of the cavity, Bq is the magnetic field strength, Ql is its loaded quality 
factor, Qa = 10^ is the 'quality factor' of the galactic halo axion signal (i.e. the ratio of their 
energy to their energy spread), pa is the density of galactic halo axions on Earth, and C is a 
mode dependent form factor given by 



C 



(18) 



where Bq{x) is the static magnetic field, E^{x)e'^^^ is the oscillating electric field and e is the 
dielectric constant. Because the axion mass is only known in order of magnitude at best, the 
cavity must be tunable and a large range of frequencies must be explored seeking a signal. The 
cavity can be tuned by moving a dielectric rod or metal post inside it. 

For a cylindrical cavity and a homogeneous longitudinal magnetic field, C = 0.69 for the 
lowest TM mode. The form factors of the other modes are much smaller. The resonant frequency 
of the lowest TM mode of a cylindrical cavity is /=115 MHz {^j^^ where R is the radius of the 

cavity. Since 10~^ eV = 27r (242 MHz), a large cylindrical cavity is convenient for searching the 
low frequency end of the range of interest. To extend the search to high frequencies without 
sacrifice in volume, one may power-combine many identical cavities which fill up the available 
volume inside a magnet's bore This method allows one to maintain C = 0(1) at high 

frequencies, albeit at the cost of increasing engineering complexity as the number of cavities 
increases. 

Axion dark matt er s earches were carried out at Brookhaven National Laboratory the 
University of Florida^^, Kyoto University , and by the ADMX cohaboration 
at Lawrence Livermore National Laboratory. Thus far, the ADMX experiment has ruled out 



KSVZ coupled axions at the nomina l halo density oi^^ pa = 7.5 10 g/cm^ over the mass 
range 1.90 < rua < 3.35 //eV^^^^^. The ADMX experiment is presently being upgraded to 
replace the HEMT (high electron mobility transistors) receivers we have used so far with SQUID 
microwave amplifiers. HEMT rece iver s have noise temperature r„ ~ 3 iC^Sl whereas ~ 
0.05 K was achieved with SQUIDs^^. In a second phase of the upgrade, the experiment will 
be equipped with a dilution refrigerator to take full advantage of the lowered electronic noise 
temperature. When both phases of the upgrade are completed, the ADMX detector will have 
sufficient sensitivity to detect axions at even a fraction of the halo density. 

The ADMX experiment is equipped with a high resolution spectrometer which allows us 
to look for narrow peaks in the spectrum of microwave photons caused by discrete flows, or 
streams, of dark matter axions in our neighborhood. In many discussions of cold dark matter 
detection it is assumed that the distribution of CDM particles in galactic halos is isothermal. 
However, there are excellent reasons to believe that a large fraction of the local density of cold 
dark matter particles is in discrete flows with definite velocities'^. Indeed, because CDM has 
very low primordial velocity dispersion and negligible interactions other than gravity, the CDM 
particles lie on a 3-dim. hypersurface in 6-dim. phase-space. This implies that the velocity 
spectrum of CDM particles at any physical location is discrete, i.e., it is the sum of distinct 
flows each with its own density and velocity. 

We searched for the peaks in the spectrum of microwave photons from axion to photon 
conversion that such discrete flows would cause in the ADMX detector. We found none and 
placed limits!^ on the density of any local flow of axions as a function of the flow velocity 
dispersion over the axion mass range 1.98 to 2.17 /ueV. Our limit on the density of discrete flows 
is approximately a factor three more severe than our limit on the total local axion dark matter 
density. 
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